Comment on 

"Dynamic Wetting by liquids of different viscosity", by 
T.D. Blake and Y.D. Shikhmurzaev 
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We comment on a recent theory of dynamic wetting, that is based directly upon a model for 
interface formation, introduced by Shikhmurzaev. We argue that the treatment of surface tension 
and its relaxation, inherent in the original model, is physically flawed. 



In a recent paper 0, Blake and Shikhmurzaev uti- 
lize a model for interface formation, proposed initially 
by Shikhmurzaev @> Q , to interpret their measurements 
of the velocity dependence of the dynamic contact angle 
at a moving contact line. They determine some of the 
parameters contained in the model, in particular, their 
viscosity dependence. In this comment, we point out that 
the key quantities upon which this model is based have 
no well-defined physical meaning. More specifically, the 
treatment of surface tension which underlies the theory 
does not correspond to any known physical mechanism. 
As a result, the model leads to consequences which are 
physically absurd and which are at odds with available 
experimental data. Furthermore, we show that the val- 
ues for model parameters, purportedly estimated in [l|, 
depart by many orders of magnitude from what is phys- 
ically reasonable. 

The model in Ref.jlJ is based on the fundamental as- 
sumption that the surface tension of a pure liquid is de- 
termined by a 'surface equation of state', which fixes its 
value as function of surface parameters (see also , sec- 
tion 2.1). As an approximation, the surface tension a is 
assumed to depend on the 'surface density' p s alone, cf. 
equations (1) and (4) of 



(1) 



where 7 and pf) are phenomenological constants. In 1], 
Shikhmurzaev's model is applied to both the free liquid- 
gas and the liquid-solid interface, but we will focus on the 
former, so there is no need for an index on a. Equation 
(TJJ makes no physical sense. Firstly, as explained in Q, 
p. 31 ff, surface density as introduced by Gibbs cannot 
be an intrinsic property of the surface of a pure liquid. It 
depends solely on the definition of what Gibbs calls the 
dividing surface between the two phases. Conventionally, 
for a pure liquid, one defines the dividing surface in such a 
way that the surface density vanishes (|^|, p. 31), a choice 
which is called the equimolar surface. Secondly, the sur- 
face is not an independent thermodynamic system that 
would allow relations between its extensive and intensive 
parameters to be defined, as is done in 0,0 El, and as 



is implied by As pointed out in H, p .33, the surface 
exists only by virtue of the bulk phases that surround it; 
it does not form an autonomous phase. Hence an expres- 
sion such as a = <j(p s ), as given in [jj, is meaningless 
as regards surface thermodynamics, no matter how p s is 
defined. 

We illustrate our criticism further by pointing out two 
of the consequences of the ansatz Q , which we think are 
absurd. Equation (13) of Q states that the surface ten- 
sion should go to zero with surface thickness h (assuming 
p s e * is constant, as is implied in the text). This is well 
known not to be the case, see p. 16 ff and p. 47, which 
provides an authoritative critique of treatments based on 
point-thermodynamics. Of course, h is not really an in- 
dependent quantity, rather, it is determined by the condi- 
tions of thermodynamic equilibrium at the surface. How- 
ever by considering, for example, different temperatures, 
different interface density profiles can effectively be re- 
alized. Given the interface profile, the surface tension 
can be calculated essentially by mechanical arguments 
0. (For a slowly- varying profile, treated in the square- 
gradient approximation of Rayleigh and van der Waals, 
the result is equation (1.43) or (3.11) of 0.) If anything, 
the surface tension will be larger for a sharp interface 
(h = 0) then for the real smooth one. The reason is 
that the real interface shape is one that minimizes the 
total free energy (grand potential) of the inhomogeneous 
fluid (see 0], p. 54 ff). Indeed, at high temperatures, 
when the interface becomes more diffuse, surface tension 
decreases (cf. Figs. 1.5, 1.6, and 6.5 0). 

Next, we consider equation (14) of 1] for the typi- 
cal time scale r, over which the surface tension reaches 
its equilibrium value after a fresh surface is created. 
The time r is claimed to be proportional to the viscos- 
ity p of the fluid, and a specific estimate is given for 
|U = 672mPas, for which the value is supposed to lie 
between r = 2.5 x lCT 6 s and 8.3 x 10~ 6 s. To a first ap- 
proximation, we in fact believe that the timescale needed 
to establish a surface tension in a pure fluid is essentially 
zero: force is transmitted with the speed of light, giving 
t ss 10 _18 s for molecular sizes. The interfacial profile, 
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and thus the surface tension, subsequently relaxes to- 
ward its equilibrium value, but even for this process the 
author's numbers appear to be a gross over- estimation. 
The thickness of the liquid-gas interface is typically 2 or 3 
molecular diameters near the triple-point (|4|, Chs 6 and 
7); this gives about 5 x 10 _10 m. Dividing this length by 
the speed of sound in a typical simple liquid, we arrive at 
t s « 10~ 12 s. In other words, r will be given by a typical 
collision time in the liquid, which is not related in any 
simple way to the fluid viscosity. 

This estimate is consistent with typical decay times of 
correlation functions in simple liquids measured by neu- 
tron scattering as well as with those found in molecular 
dynamics simulations ||. Note also that Brillouin scat- 
tering from density fluctuations in bulk liquids, where 
the wavelengths involved are much larger, of the order 
of 5 x 10~ 7 m, is characterized by shifts of frequency of 
typically \Q w Hz and that Brillouin spectra for liquid in- 
terfaces correspond to the same range of frequency shifts 
0. We are not aware of any physical mechanism that 
would give rise to relaxation times in the order of 10 _6 s. 

One might argue that as the authors' estimates are 
based on their own measurements of a moving contact 
line, they should provide independent evidence for the 
consistency of the assumptions of the model. Unfortu- 
nately, we do not believe this to be the case. The reason 
is that the authors' theory refers to what they term the 
dynamic contact angle 9d, determined solely by a balance 
of surface tensions, cf. equation (7) of 1]. Viscous forces 



and, therefore, interface bending does not enter their de- 
scription. The measurement of the interface angle was 
however performed at a scale of about a mm. It is well 
known that the interface near a moving contact line is 
highly curved 0], which is the result of viscous forces 
which therefore cannot be ignored. This is best appre- 
ciated in the case of a perfectly wetting fluid, where the 
contact line is preceded by a precursor film jjj . Hence no 
interface formation is taking place, yet on macroscopic 
scales measured contact angles have a speed dependence 
consistent with Tanner's law |7J . This implies that in the 
partially wetting case, considered here, any effects of in- 
terface bending would have to be carefully subtracted for 
a correct interpretation of experimental data. 
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